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Abstract 

We calculate the vacuum expectation value of the stress-energy ten- 
sor for a quantised bulk scalar field in the Randall- Sundrum model, and 
discuss the consequences of its local behaviour for the self-consistency 
of the model. We find that, in general, the stress-energy tensor di- 
verges in the vicinity of the branes. Our main conclusion is that the 
stress-energy tensor is sufficiently complicated that it has implications 
for the effective potential, or radion stabilisation, methods that have 
so far been used. 



1 Introduction 

It is widely being considered whether our spacetime has more than four 
dimensions. Higher dimensional models have arisen in the context of the 
search for a unified field theory since the proposals of Nordstrom, Kaluza 
and Klein ^ in the early years of general relativity. These ideas have been 
extended over time, and have recently been revived |2j in the form of the 
so-called "brane- world" models, in which the (3-1-1 )-dimensional universe we 
observe exists as a membrane, or 3-brane, in a higher dimensional spacetime. 
In the framework of string theory, Horava and Witten jBj have constructed 
such a model that may accommodate a realistic gauge theory. 
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In 1999, Randall and Sundrum (RS) [1] proposed that a scenario closely 
related to the Horava-Witten model could naturally explain the problem 
of the large hierarchy between the fundamental scale of quantum gravity, 
the Planck scale Mpi ~ 10^^ TeV, and the scale of the standard model 
of particle physics, the electroweak scale Mew ~ 1 TeV. In this scenario, 
two flat, parallel 3-branes of opposite tension sit at the fixed points of a 
/7j2 orbifold. The "bulk" between the two branes forms a slice of 5- 
dimensional anti-de Sitter spacetime. The geometry is then non-factorisable, 
the four-dimensional part of the metric being multiplied by a "warp" factor, 
an exponential function of the fifth dimension: 

ds'^ = e-'^"^y^r]^^dx^'dx'' - dy^, (1) 

where = — {k ~ Mpi) and — vrr < y < vrr is the extra dimension. 
The branes are located at y = and y = nr. The brane at y = vrr known 
as the visible brane and represents the universe as we know it. Conversely, 
that at y = is called the hidden brane. In the original RS model, gravity is 
the only field present in all five dimensions. All of the standard model fields 
are confined to the visible brane. Alternatives to this somewhat artificial 
restriction have been investigated by many authors (see e.g. 0)- It is the 
possibility of bulk matter fields, in particular scalar fields, that we deal with 
in this paper. 

The remarkable point about the RS model is that any mass parame- 
ter mo on the visible brane in the fundamental higher dimensional theory 
will correspond to a physical mass of e~'^'^^mQ when measured in four di- 
mensions. By taking A;r ~ 12, it is possible to generate a TeV mass scale 
from the Planck scale, without introducing any large hierarchies between 
the fundamental parameters of the theory. 

However, in the basic model, the compactification radius r is left unde- 
termined. In order to make the model physically acceptable, a mechanism 
for stabilising r must be found. A classical mechanism has been suggested 
in 0. 

In the context of the older Kaluza-Klein theories, Candelas and Weinberg 
[Jj in 1984 found that quantum effects from matter and gravity fields could 
fix the size of the extra dimensions in a natural way. Various authors have 
investigated whether a similar procedure could fix the value of r in the RS 

model iBunmni. 

In a detailed analysis was given of the requirement of the self- 

consistency of the model, that the model be a solution of the quantum- 
corrected Einstein equations, showing that the problem is not just one of 
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minimising the vacuum energy density of the bulk field. The authors found 
that it is possible to find self-consistent solutions that solve the hierarchy 
problem only at the price of severely fine-tuning the brane tensions. One 
of the present authors is currently applying their method to the case of 
de-Sitter branes, which is of interest to inflationary cosmology |12j . 

These, and indeed most, calculations have been performed by introduc- 
ing a bulk scalar field to the model and reducing the theory to four dimen- 
sions by summing or integrating over the extra dimension. However, impor- 
tant effects can be overlooked in this approach. In this paper, we investigate 
the local behaviour of the full stress-energy tensor for a quantised bulk scalar 
field, and discuss what the consequences are for the self-consistency of the RS 
model from this standpoint. We show that the problem of self-consistency 
and stabilisation of the RS model is more complicated than has generally 
been assumed. Other related work includes ITU [T^ . 

2 Set-up and general method 

For ease of calculation, we first rewrite the metric in Riemannian form 
so that all momenta will be positive-definite: 



where a caret denotes a (D + l)-dimensional quantity: fi = 0,1,. ..,D, so 
that = y, whereas = 0, 1, [D — 1). 

Into the bulk of the above background spacetime, we place a single real 
scalar field (p of mass m. We allow a general coupling to the scalar curvature 
R with dimensionless constant The action S of the field is 





X I det g 



(3) 



The stress-energy tensor for such a scalar field is given by 




(4) 



where G^p is the Einstein tensor, which is given by 



(5) 
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away from the branes. 

Rather than working directly with with this expression for the stress 
tensor, we re-express the stress tensor as a vacuum expectation value in 
terms of the Green function G{x, x') for the field. 

Expressed as a D-dimensional Fourier transform, the Green function for 
a scalar field in a slice of {D + l)-dimensional anti-de Sitter spacetime is 



G{x,x') 



(2^ 



(6) 



with 



Gp{y,y') 



vr D 
— e 2 
2k 



((T>+(T<) 



jy{ip/ak)Hl{ipe''> /k) - hl{ip/ak)Jy{ipe''> /k) 



jy{ip/ak)hl{ip/k) - hl{ip/ ak)jy{ip/k) 
'j,{ip/k)Hl{ipe^</k) - hl{ip/k)Mipe''</k) 



(7) 



where is a Bessel function and Hi, is the Hankel function of the first kind, 
both of order v. We have defined: 

ju{z) = ^D{l-AOMz) + z^Mz), 

Kiz) = \d{1-A0HI{z)+z^HI{z), 



and have written a 



and 



Now consider 



-knr 



0-< 



We define 

_ f o-(y) for y > y' 
" 1 o-(y') for y <y' 



a{y') for y > y' 
a{y) for y < y' 



G{x,x') = (O|0(£)0(x') 
Using this, together with it can be shown that 

= G'(X,X')U'=£, 



(8) 

(9) 

(10) 
(11) 



4 



{d^(l)d^(/)) 

{dy(t)dy(t)) 

{df,(pdy(p) 
\yy) 



{{^%,y) = 
We also note that 

d^d,yG{x,x')\s:r 



0, 

-ke~'^''5^u{dy + dyi)G{x,x')y=x, 
{2dydy> + dy + 9y)G'(x,x )U'=x, 
0. 



p'^Gp{y,y')\y>-_ 



(12) 



(13) 



The vacuum expectation value of the stress tensor @ can therefore be writ- 
ten in terms of the Green function as 



flU/ 



(2^ 



1 



1 



D 



D 



2 



^{D-l)ke-^^{dy+dy' 



Gp{y,y')\ 



y =y^ 



yyi 



(2^ 



ldydy,-^e^-p^+^j{D-l)k' 



1 



— ^Dk{dy + dy: 



Gpiy,y') 



\y'=yi 



(14) 



where we have used equation ^ for the Einstein tensor. The off-diagonal 
elements are zero. 

Substituting Gp from Q into the above expression gives the stress tensor 
for a bulk scalar field in the Randall-Sundrum model before renormalisation. 
Care must be taken in taking the derivatives of Gp{y,y'). We use the sym- 
metrised relation 



^Gp{y,y')\ 



y =y 



hm -nGp{y ■> v')\y'<y + Gp{y^ y')\ 



y'^y 



y >y)^ 



(15) 



where & is a differential operator. 

The expressions for the stress tensor components so obtained are very 
complicated in general. To perform the integrals over p and to carry out the 
necessary renormalisations, numerical methods must be used in all but the 
simplest case of massless, conformally coupled fields. 
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3 Massless, conformally coupled fields 

The case of massless, conformally coupled fields is especially simple. The 
Bessel functions are of order v = 1/2, so that they become elementary 
trigonometric functions, allowing the stress tensor to be calculated exactly. 
This provides a useful check on the method. 

For the conformally coupled case, m = and ^ = (D — 1)/4Z), so that 
we obtain: 



Gp{y,y') 



gii^(<x>+a<) 

psinh ^{e^'"'^ — 1) 



P 



P, 



cosh -r{e'^'^ — 1) cosh — (e' 



kTrr\ 



k 



k 



(16) 



for the Fourier transform of the Green function. This is the same as the 

(-D-1) / , ^ 

corresponding function in flat space multiplied by e 2 l'^>+o'<i. 
The stress tensor ()14() takes the form 



1 1 
D~2 



+■ 



D~l 



yyi 



AD 
4D 

(2^ 
{D-1 



ke-^^idy+dy') 



Gp{y,y')\y'=y , 

k' 



■k{dy + dyl) 



Gp(y,y')l 



y =y • 



(17) 



Substituting (fTO)) into (fT7|) . we find that 

1 



^pcoth^(e'=- 



1 



yyi 



(2^)^ " k 
(0,.coth|(e--l). 



1), 



(18) 



The stress tensor is therefore traceless, as it should be in the conformal case, 
so we only need to consider the (0,0) component. First, we switch to polar 
coordinates and write 



00/ 



1 

'2D^ 



1 



2D-It,d/2y^D/2)' 

r dpp^ coth^(e'=^''-l). 
JO k 



(19) 



6 



The integral / requires regularisation, as coth ^{e^'^^ — 1) — > 1 as p — > oo 
First, we rewrite the integral as 

2g-2/3p/fc \ 



CO 



1=1 ^pp^ 1 1 + J , (20) 
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where we have put j3 = — The first term in the integrand is responsible 
for the divergence. We employ a y-dependent cutoff e~^A, to obtain 

= D + l + &(« + !), (21) 

where Ci? is the Riemann zeta function (see jl7j^. 

By substituting the first term into (|19|) . we can see that the diver- 
gent part of (Too) is proportional to e~'^^A. This can be absorbed into 
the (0,0) component of the Einstein tensor 0, and is the reason for in- 
cluding y-dependence in the cutoff. Intuitively, this comes from the en- 
ergy/momentum scaling due to the warping in the extra dimension. In 
making the i^-dimensional Fourier transform (jH)), we have singled out the 
extra dimension. The coordinate y measures coordinate length in this di- 
mension, rather than proper length, the difference involving the warp factor. 
Something similar was performed in JOl) in which a rescaled cutoff was used 
on the y = nr brane, from a D-dimensional point of view. 

After these considerations, we can drop the divergent term and obtain 

{f00)ren = - ^VT" (4.)-^/^ (l^)''^' ^^^^^^^^ ^^^^ + 

(22) 



Similarly, 



{fyy)ren = f ^"^/^ (4vr)-^/2 (l^) ""^^ ^^""^'^^^ (^^) ^^^^ + 

(23) 

This result, for the stress tensor for a massless, conformally coupled 
scalar field, agrees with what would be obtained by performing the calcu- 
lation via a conformal transformation to flat spacetime. Bearing in mind 
that limfc_>o = and setting D = 3, we obtain the well-known scalar 
Casimir result for parallel plates in a (3 + l)-dimensional flat background: 

^2 



where we have put i = irr. 
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4 General mass and coupling 



As mentioned above, the calculation can in general only be carried out 
numerically. To begin with, we choose particular values of interest for the 
dimension, mass and ^, and then obtain the forms of the Green function Q 
and its derivatives for those values. These functions are then substituted 
into the integral expression for the stress tensor 1)14(1 . We also set A: = 1, so 
that we work in units of the warp scale, which is of the order of the Planck 
scale. 

As before, the integral (|14|) for {T^(,) is ultra-violet divergent. In order 
to renormalise the integral, the asymptotic behaviour of the integrand for 
large p needs to be found. 

First, we make the replacements ^H] 

Ju{iz) = i^Iuiz), 

Hl[iz) = -^K,{z), (25) 

where ly and Ky are the modified Bessel functions, so that we can use the 
following for large z jllj 

Iy{z) = -^S^(^), 
y/znz 



with 



where 



and 



Ky{z) = ^/_e-S^(z), (26) 



S^(z)=.^afcZ-\ (27) 

fc=0 

(_i)fcr(.. + fc+i) 
"'^ = 2^fc!r(.-fe + i)' ^''^ 



S-^(z) ~ S^(-z). (29) 



In practice, we take the first six terms in the sum ()27j) . 

At this stage, we choose a value for a, and use the Maple VI asympt 
function to find the large p behaviour. In this way, we obtain the asymptotic 
behaviour of the integrand in (|14|) . for which we write t(p), as a power series: 

, d\ do d-^ di 

t p ~ap + (io + — + ^ + ^ + 4 + - (30) 
p p^ p-^ p* 
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Now, we have to calculate an integral of the form 



d^pt{p) 



dpp^~^t{p). 



(31) 



Therefore, in order to render the integral finite, we need to subtract off terms 
of the asymptotic series (|30|) up to and including the term of order 1/p^ . 

However, this introduces a problem: the final term in the subtraction 
diverges for small p. This is remedied by multiplying this term by e~'^/^, 
where ^ will be the renormalisation scale. 

The value of the integral is then dependent upon ^i. We now add and 
subtract a term {d£)/p^)e~^^'^ . This allows us to isolate the renormalisation 
scale dependence: 



dpp^~H{p) 



dpp^ ^ [t{p) — ap — do 



di 
P 



where we have used 



roo 1 

/ dp-{e~^/P -e-'^/P) = \n 
Jo P 



+ dDlnfi, (32) 



(33) 



We then calculate the renormalised integral (|32|) numerically for a set of 
y- values. 

In order to test the numerical technique, we have used it to re-calculate 
the stress tensor in the conformal case. We have chosen a = 10~^, as 
although a "realistic" value would be a ~ 10~^^, a smaller negative order 
of magnitude allows the computer to calculate effectively and allows the 
features of the results to be seen more readily, without qualitatively changing 
the results. This value of a places the branes at y = and y = 2.302585.... 
We also take D = A, m = and ^ = 3/16. 

The results are plotted in Figure ^ together with the exact analytical 
expression H23|) . It can be seen that the numerical results agree perfectly 
with the exact calculation. 

Having established that the method works, we move to the massless, 
minimally-coupled case. This is the next simplest case, though already the 
calculation cannot be performed exactly. We take the same parameters as 
above for a, D and m, though we now set ^ = 0. 

The results are plotted in Figure [21 It is evident that the stress tensor 
diverges near the branes. We have also examined other cases with various 
values of m and ^, obtaining similar results. 
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1 1.5 

y 



Figure 1: The conformally-coupled case in (4+1) dimensions. The perfect agreement 
between the exact expression for {Tyy)ren (dashed Une) and the numerical data (points) 
is illustrated. The branes are located at y = and y ~ 2.3. fc is set to 1 and a — 0.1. 



4.1 Surface divergence 

In |19j . Deutsch and Candelas have argued that surface divergences in the 
Casimir effect arise from the high wavenumber (momentum) behaviour of 
the Fourier representation of (7/ii>)- Taking this hint, we use it to help us 
isolate the part of the integrand t{p) responsible for our brane divergence. 

We again take our expression for t{p) in terms of the functions (|26() . 
This time, we only take the first term in the sum (|27|) . We obtain a sum 
of terms of the form exp(/(p, y)), where n is some power and f{p,y) is 
some function in which p is linear. We pick out the terms with the highest 
power n. Then, we notice that f{p,y) falls into two categories concerning 
its y-dependence: either it is linear in y or contains e"^*-^^. 

We need to consider what happens to each of these kinds of term when 
integrated over momentum space. The first simply give either terms diver- 
gent for all y, part of what is removed by renormalisation, or give terms 
finite for all y. The second kind of term, however, is more complicated. 
These terms will turn out to be responsible for the brane divergence. We 
therefore discard the terms of the first kind. 

After collecting the remaining terms and simplifying, it turns out that 
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we are left with just two terms: 



tipUv ~ e^^-^^> ^e^K^^-Vaj + g^PU-e^ij . (34) 

The first term is responsible for the y = irr brane divergence and the second 
for the y = brane divergence. Integrating over p space, and taking the limit 
as y approaches the appropriate brane, we find the following expressions for 
the near-brane behaviour: 

- a 

{Too)ren\y^O = '^p+l^ 



where 



a 



.l)D+i2-(D+2)^-(D+i)/2 _ _ ^) p ('^^ ) . (36) 



The dashed lines in Figure |21 represent these expressions in the massless, 
minimally coupled case, which fit the numerical results well close to the 
branes. Figure 01 shows the fractional deviations of the numerical data from 
the near-brane expressions. 

It should be noted that the divergences only disappear for ^ = i.e. 
in the conformally coupled case, consistently with what we found earlier. 

At face value, our result that the bulk stress-energy tensor diverges non- 
integrably as one approaches the branes would seem to contradict the finite 
result for the total vacuum energy as calculated in in which the renor- 
malised sum of the zero point energies for each mode of the bulk field was 
calculated. 

Similar problems regarding the Casimir effect between conducting sur- 
faces in flat space were first analysed in detail by Deutsch and Candelas 
They argued that this mode-sum energy is generically different from 
the integral of the renormalised energy density (Too)ren) because the op- 
eration of forming the volume integral does not in general commute with 
that of renormalisation. Following this analysis, Kennedy, Critchley and 
Dowker |2nj showed that the resolution lies in renormalisation of surface 
(on-boundary) terms in the effective action, which does not allow a local 
interpretation. This work has recently been extended to generally-coupled 
scalar fields for general bulk and boundary geometries |21j . 

Essentially, our stress-energy tensor is strictly-speaking incomplete: the 
full renormalised tensor is identical to ours, but contains an infinite term on 
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-6 I 1 1 ' 1 1 

0.5 1 1.5 2 2.5 

y 

Figure 2: The (4+l)-dimensional massless, minimally-coupled case. The points show 
the data for (Too)reri from the numerical calculation. The dashed lines represent the 
analytic expressions for the leading near-brane behaviour. The parameters a and k and 
the positions of the branes are the same as for Figure 

the boundary. On integration, this surface term exactly balances the infinite 
contribution from the local divergences in the bulk part of the stress-energy 
tensor. The result is then equal to the sum of the zero-point energies. 

5 Discussion and conclusions 

In this paper, we have calculated the stress-energy tensor for a quantised 
bulk scalar field coupled non-minimally to the curvature in the background 
spacetime of the RS model. While an analytic result is not possible in 
general, the case of a massless, conformally-coupled field is calculated ex- 
actly. In the general case, we have produced numerical results that reveal 
the stress-energy to diverge near the branes, and have obtained analytic 
expressions for the leading behaviour of these divergences. 

Self-Consistency 

Self-consistency requires that the RS metric ((21) be a solution of the quantum 
corrected Einstein equations. This means including the bulk field stress- 
energy tensor on the right hand side of the Einstein equations with the 
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Figure 3: The deviations from the near-brane behaviour. The difference between the 
near-brane curves and data points in Figure|21are plotted as a fraction of the data, showing 
the convergence as the branes are approached. 

RS metric as an ansatz (with general cr{y)) and solving for <T(y). For the 
original RS metric to be strictly self-consistent, the result for a{y) should 
be of the form k\y\. The constant k might not take its classical value, but 
may include a small quantum correction, this correction being the same for 
every component of the Einstein equations. 

In the case of conformal coupling, it can be observed immediately that 
the traceless stress-energy tensor has a different form than the Einstein 
tensor, and hence each component will give a different correction to k, if 
indeed a RS-type solution can be found. This in itself implies a lack of 
strict self-consistency. However, if the stress-energy tensor is small (this is 
in fact the case: e.g. |(7bo)ren| ~ W~^^k^ at most for D = 4), then we might 
expect the RS solution to be approximately self-consistent, in the sense that 
(T(y) = k\y\ plus extra small terms due to quantum effects, which may differ 
for each component. The true self-consistent solution might then not differ 
from the original RS model significantly. 

In the general case, the stress-energy tensor diverges close to the branes. 
This would also be expected in the conformal case if the branes are curved, 
however slightly . Although the generic stress-energy tensor is not trace- 
less, it would also not be expected to have to same form as the Einstein ten- 
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sor, except possibly in some special case. Therefore, any argument relying 
upon smallness of the stress-energy tensor will not apply, and a conclusion 
is then that, generically, the RS metric is not self-consistent, even approx- 
imately, in the presence of a bulk scalar field. If there is a self-consistent 
solution, then a{y) will not approximate to k\y\. 

As mentioned above, the Casimir effect also suffers from divergences. 
However, the imperfect nature of the conducting surfaces becomes apparent 
(in the form of a finite skin-depth |19|) and destroys the divergence long be- 
fore the backreaction of the fields between the plates becomes measurable. 
In our case, the analogue is the quantum "imperfection" of the branes them- 
selves, on the order of the Planck scale, at which quantum-gravity effects 
would come into play and alter the situation. Due to the 1 /y^^^ behaviour 
of the divergences, the backreaction should become significant before this 
region is reached. 

It should be noted that this conclusion might not hold in a supersym- 
metric theory, where a fermionic contribution could exactly cancel the bulk 
field vacuum energies. We have also not taken into account the effect of the 
surface terms in the stress-energy tensor. Furthermore, it would be inter- 
esting to examine the consequences of introducing a potential to the scalar 
field action. 

Whether or not the RS model with bulk fields can be rendered self- 
consistent, the main conclusion of this paper is that the stress-energy tensor 
is sufficiently complicated that it has implications for the effective potential 
methods that have so far been used. The important difference with the older 
Kaluza-Klein scenario of Candelas and Weinberg [7] is that assuming spheres 
for the extra dimensions, as in that case, results in a homogeneous space 
for which the stress tensor is determined by symmetry from the vacuum 
energy. For the RS model, homogeneity is lacking so that it is not necessarily 
sufficient to use an effective potential in general. 
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